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INVESTIGATION OF THE LIFT DISTRIBUTION 
OVER THE SEPARATE fflNGS OF A BIPLANE* 
By D. Kuchemann 



An investigation is made of the mutual interference 
of the wings of a "biplane under the general assumption 
that each wing may be replaced by a vortex system of the 
type given by the Prandtl wing theory. The additional ve- 
locities induced at each wing by the presence of the other 
are determined by the Biot-Savart law and converted into 
an equivalent change in angle of attack, the effect being 
that of an additional twist given to the wings in changing 
their lift distributions. The lift distributions computed 
in this manner for several airplane types are compared 
with the results of measurement. 



Notation 

In general the subscript o will denote the value 
for the upper wing, whereas the subscript u will denote 
the corresponding value for the lower wing; 

y, z, rectangular coordinates (right-hand system); the 
y axis laid in the direction of the bound 
or lifting vortex of the lower wing positive 
to starboard, origin taken at the center of 
the wing, and the x axis forward 

a, stagger of the biplane system = separation of 
lifting vortices in the x direction 

h, wing gap = separation of lifting vortices in the 
z direction 



t , wing chord 
b , wing span 



* "Berechnung der Auf triebsverteilung uber die einzelnen 
Fl-iigel eines Doppeldecker s . " Luf tf ahrtf or schung, Oct. 10, 
1938 , pp, 543 - 55 . 
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T| = 2y/b, nondimensional coordinate in the y direction 

a and a , angle of attack measured from the zero 
lift direction of the upper wing 



cci , self-induced angle of attack 



x ou 



angle of attack induced at the upper wing by the 
presence of the lower wing 



a iuo 



correspondingly for the lower wing 
a e . effective angle of attack 



f , effective rise of camber 
e ' 

V, undisturbed velocity at' infinity 

u, w, induced velocity components in the x and y direc- 
tions, respectively 

A, W, M, lift, drag, and moment per unit length 



c a 
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section lift coefficient 



section drag coefficient 



c = — B — , section moment coefficient 
m 



c = , total lift coefficient of a wing 

a, " 



c a d> 



total lift, coefficient of cellule referred to the 
total area. 



P, circulation 
r c a * 

7 = = , nondimensional coefficient of the circu- 

b V 2 b 

lation 

2 



b /? , aspect ratio of a wing 
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I.' REVIEW OF PREVIOUS.. WORK DONE OH THE BIPLANE PROBLEM 



The solution of the monoplane wing problem' having been 
obtained, that of the biplane lifting system consisting of 
the combination of two. airfoils presents new problems for 
the airfoil theory to solve, because of the fact that each 
wing is located in the disturbed flow produced by the other. 
Each wing of the cellule therefore exerts an essential in- 
fluence on the position and magnitude of the aerodynamic 
forces on the other. From the practical viewpoint, an an- 
swer is required to the following two questions: 



1. How large is the total lift of each of the wings 

and of the biplane combination? 

2. How is the lift distributed over the span of each 

wing? 

In connection with the second question, a third one 
arises, namely, 

3. What lift distribution gives the minimum induced 

drag of the lifting system? 

The first question has already received the attention 
of a large number of investigators, whereas an attempt to 
answer the second will be the object of the present paper. 
To the third questions a general solution had already been 
found very early in the development of the theory. 

I should like first to give a brief review of the 
work that has so far been done on the biplane problem. 
This work is almost entirely concerned with the answer to 
the first question. It is convenient to distinguish two 
groups of theoretical investigations. A first group con- 
siders wings of finite span (three-dimensional problem). 
In those investigations in which only the total lift is to 
be determined, assumptions are made on the' lift distribu- 
tion so that each wing is replaced either by a horseshoe 
vortex or by a vortex system with elliptic lift distribu- 
tion. It is then possible to consider theoretically the 
mutual interference of such vortex systems. The finite 
chord of the airfoil is only approximately taken into con- 
sideration. In contrast to this group there is a second 
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larger group of investigations, which takes the effect of 
the finite chord into account by considering the two di- 
mensional flow, "but which does not take into account the 
effect of the finite span or only estimates it. This sec- 
ond mode of treatment has the advantage over the first in 
that the mathematical aids of conformal transformation and 
integral equations are to a large extent made available 
and this also explains the large number of these investi- 
gations. The method has the disadvantage that the effect 
of the chord, which is mainly taken into account, is in 
general considerably smaller than the effect of the finite 
span of the wings. We shall first consider the group of 
inv e st igat i ons , whi ch deals the solution of the first prob- 
lem, namely, that which takes the three-dimensional prob- 
lem as the starting point. 

The first computation of the mutual, interpf erence ef- 
fect of two airfoils which led to practical formulas is 
due to A. Betz (references 1 and 2). The general method 
of treatment is essentially the same as that given in the 
present paper and will be presented in detail, below. As 
a basis, the assumption was made that the lift is uniform- 
ly distributed over the span (simple horseshoe vortex). 
It was possible in this way to determine the change in 
total lift to a first approximation. A refinement of the 
method by L. Prandtl (references 3 and 4) assumes elliptic 
distribution for both wings. • The • carrying through of the 
computation was in this case, however, made considerably 
more difficult so that only simpler problems than those 
considered by A. Betz could be solved. In connection with 
this, H. von Sanden (reference 5) gave a numerical compu- 
taion of the deviation in the obtained mean downwash veloc- 
ity induced by the trailing vortices as given by the theory 
when uniform and elliptic lift distributions, respectively, 
were assumed. There was obtained the oft-mentioned result 
that the differences in those mean values that arc of im- 
portance for the total lift and drag are not usually very 
large. An attempt was made later to take greater account 
of the wing chord. The due consideration of the wing chord 
is, as we shall see later, important for the reason that 
not only is the velocity at one wing changed by the pres- 
ence of the other but also the streamlines in the neigh- 
borhood receive a curvature. H. K. Bose (reference 6) has 
particularly investigated the effect of curvature of the 
streamlines. A comparison of his results with those ob- 
tained for the streamline curvature in the present paper 
showed good agreement. A comprehensive presentation of 
this method of treating the biplane problem has been found 
by E. Fuchs (reference 7), who has developed and refined 
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the method in a paper as yet unpublished:, so that the total 
lift may be determined by this method with sufficient ac- 
curacy inmost of the practical cases occurring. This es- 
sentially completes the list of the first group of papers, 
which concerns itself directly with the three-dimensional 
problem. All these theories are based on particular as- 
sumptions on the lift distributions over the two wings and 
compute the change in the total lift. 

A brief reference may be made here to the papers which 
deal with the answer to the third question and which are 
also concerned with the three-dimensional flow about the 
biplane cellule. These are the papers by L. Prandtl (ref- 
erences 3 and 4) and by M. M. Munk (-reference 8) in which 
Is determined the lift distribution which gives the. .mini- 
mum induced drag of- the biplane cellule. These papers also 
contain general theorems on the induced drag of biplane 
lifting systems. 

Again returning to the reports that have appeared on 
the solution of the first of the above-mentioned problems, 
we come to that group in which the problem is reduced to a 
two-dimensional one. It is not possible here to give a 
detailed account of all the methods used. Essentially the 
two-dimensional flow about two infinite airfoils was inves- 
tigated with the aid of the theory of conformal transforma- 
tion or that of integral- equations. The first investiga- 
tion of this kind with reference to practical application 
is that of ff. Kutta (reference 9), whose investigation 
was again carried out more in detail by H. Grammel (refer- 
ence 10). In the papers by E. Pohlhausen (reference 11), 
M. M. Munk (reference 12), 33 . Eck (reference 13), T. Moriya 
(reference 14), H. Glauert (reference 15), M. Lagally (ref- 
erence 16), C. B. Millikan (reference 17), I. Tani (refer- 
ence 18), C. Ferrari (reference 19), E. Pistolesi (refer- 
ence 20), G. Schmitz (reference 21), P. Teofilato (refer- 
ence 22), J. Bonder (reference 23), P. Koebe. (reference 24), 
I. E. Garrick (reference 25), and E. Graeser (26), there is 
treated principally the conformal transformation of two cir- 
cular arcs into profiles, and conclusions are then derived 
for the biplane. The amount of computation for specified 
profile shapes is, however,, quite considerable. In some 
of the papers the very important effect of the trailing 
vortices (see computation under section III) is sufficient- 
ly well taken into account,' B. Eck (reference 13) replaces, 
for example, the effect of the' trailing vortices by a uni- 
form downwash velocity and thus the flow about the biplane 
by a two-dimensional flow, a mean circulation, and an 
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elliptic lift di stribution being assumed for "both, wings. 
In various French papers the results of the two dimension- 
al theory are given and partly extended to the three-di- 
mensional ca^es. A comprehensive report covering this 
work by Li. Nenadovitch will be found in reference 27. He 
also conducted a large number of measurements on approxi- 
mately two-dimensional flows. This practically completes 
the list of theoretical investigations. 

The theoretical papers are supplemented by a number 
of exper imenta-1 papers. In addition to the above-mentioned 
French measurements (reference 27), there are those of M. M. 
,4unk (references 23 and 37); also, a large number of English 
measurements of which we may mention those of W. L. Cowley 
and H. Levy (reference 29), W. L. Cowley and C. H . H. Lock 
(reference 30), and W. L. Cowley and L. Jones (reference 31) 
(some further though not so detailed measurements are found 
in the British A.S.C. Reports and Memoranda); further, the 
results of the Aerodynamic Experimental Institute of Got- 
tingen (reference 32) and those of F. H. Norton (reference 
33), a quite comprehensive series of measurements by the 
Toronto University, Canada (references 34, 35, and 36) in 
which are included also pressure distribution measurements; 
and finally the tests of A. I. Fairbanks (reference 38), 
0. E. Loeser (reference 39), tf. Knight and E. W. Noyes 
(•references 40 and 41), and S. A. Pearson (reference 42), 
Some of these measurements will be used for checking our 
results. Unfortunately, measurements of a kind that will 
satisfy all requirements and from which the lift distribu- 
tions over the wing spans may be obtained exist only in 
very small number, so that a systematic test program includ- 
ing pressure distribution measurements in a large wind tun- 
nel would be very desirable. 

In conclusions reference may be made to a group of 
papers which are concerned with the development of working 
charts for practical application, utilizing the theoretical 
and experimental results available. Of chief importance 
are those of .P. Kuhn (reference 43) and W. S. Diehl (refer- 
ence 45). In these a knowledge of the total lift of the 
cellule is assumed and this lift is then distributed be- 
tween the upper and lower wing. A determination of this 
kind can be only very rough since the parameters of stagger, 
gap, span, chord, lift slope, angle of attack, and decalage, 
as we shall see below, do not always enter linearly. A 
comparison of our r e suit s. wi th those of W. S. Di ehl , (refer- 
ences 44 and 45) indicated tne possibility of making the re- 
sults of the latter agree accurately with ours. 
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II. STATEMENT 03* THE -PROBLEM 



. In the present report an attempt will be made' to an- " 
swer the second of the above-mentioned questions (section 
I). We shall therefore consider a "biplane of finite span 
and try to compute the lift distribution over each wing 
for a given shape of the biplane cellule. The exaqt' dis- 
tribution is important for the static computation of the 
cellule and for the . computation of the field of flow about 
the wings, particularly for the dpwnwash at "the tail; and 
finally, for the determination of those pos.itibns along 
the span at which wi th ' incr ea sing angle of attack the' max- 
imum/lift coefficient c amo ^. is first attained. 

For the main portion of the computation we shall re- 
place each of the airfoils by a straight lifting vortex 
filament with variable circulation as is done in the 
Prahdtl theory for the monoplane (fig. 1). The lifting ■ 
line will be passed through the center of the profile sec- 
tion. It will, however, appear below that, particularly 
for small gaps and in the range near the zero angle of at- 
tack, this simple assumption will not suffice and it will 
be necessary to replace each wing by at least two vortex 
filaments. How this is to be done will be shown more in 
detail later, section V. This substitution takes account 
of the finite chord of the wing. We also take into consid 
eration the wing chord in all cases in the determination 
of the effect of the streamline curvature. 

We shall first determine the mutual interference of 
the two wings of the biplane. The vortex system of the 
lower wing induces' at the location of the upper wing addi- 
tional velocities which may be converted into a change in 
the angle pf attack of the upper wing, and conversely. 
This is equivalent to a twist given to each' wing by the 
presence of' the other. With the aid of this additional 
angle of attack, which we shall denote as "externally in- 
duced," and the ' self -induced angle of attack,' we can then 
compute the effective angle of attack from the geometrical 

cce 0 = <Xo - cc io - a i<m (1) 
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The subscripts o and u : .will denote respectively the up- 
per and lower wings, so that, for example, 0Ci ou is the 

angle of attack of the upper wing induced by the presence 
of the lower. With the aid of these equations, the two 
problems set up by L. Prandtl for the monoplane can be 
solved for the case of the biplane. The first problem; 
namely, to determine the wing shape for a given lift dis- 
tribution, offers no difficulties. In this paper we shall 
consider the practically more important problem, ■ namely , 
for a given biplane cellule to compute the lift distribu- 
tion over each wing. We must first determine the external- 
ly induced or interference angles of attack a* and 

ou 

CCi as functions of the shape of the biplane cellule' and 

"HQ 

the angles of attack of each wing... 

III. DETERMINATION OF THE INTERFERENCE ANGLES OF ATTACK 
AND THE SETTING UP OF -THE GENERAL EQUATIONS 
FOR THE CIRCULATIONS 



The interference angles of attack cu and a,.; 

■'■ou x uo 

result from the combined effect of many factors, which will 
now be considered. In general, we determine according to 
the Biot-Savart law the velocities induced at the midpoints 
of the wing sections of one wing by the vortex elements of 
the other and convert the x and z components u and w, 
respectively, in the usual manner into an increment in the 
angle of attack. 

The bound or lifting vortex of each wing induces at 
the location of the other wing a change in angle of attack: 

1. Through a change in the undisturbed longitudinal 

velocity ( x component) 

2. Through an additional upwash or downwash component 

(z component) 

3. Through curvature of the streamlines as a result 

of the change in w with x. 

The trailing vortices induce at the wing a change in 
angle of attack: 
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1. Through additional upwash or downwash velocities 

(z component ) 

2. Through curvature in the streamlines due to change 

of w with x. 

We first consider the upper wing (fig. 2) . The de- 
tailed computation of the velocity components according to 
the Biot-Savart law may "be dispensed with here (see for 
example, Glauert (reference 15)). We need only consider 
somewhat more closely the conversion into angle of attack, 
the effect of the change in the undisturbed longitudinal 
velocity, and the curvature effect. 



In order to determine the first, we start with the 
Prandtl airfoil equation for the monoplane ( * which may, for 
example, be the upper wing alone: 

b. 



+ -2 



lo i p 2 a r 0 d y 0 « 

= a 0 / - (3) 

▼o*o > 4 tt V Q ,L i y Q > y Q - y Q « 



and assume that 



V = V + A V = V (1 + 



where V is the undisturbed velocity and A V its incre- 
ment, that is, the x component of the velocity induced 
by the lifting vortex of the lower wing at the upper wing. 
Equation (3) thus becomes oh multiplication by (1 + AV/v) 

b, 

+ 

r 

o AV 1 



a 0 + cc f 




.Cj Tt fl " "V 4 tt V 
o ° 



The velocity increment ■ A V thus corresponds to a change 
in the angle of attack 



*I't is assumed that the flow is potential. 
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A V 

A a 0 = cc 0 — 

( see f ig. 4) . 

The curvature of the streamlines due to the lifting 
vortex of the lower, wing gives rise to an angle of attack 
increment (see, for example, Giauert , reference 15) 

A cc 0 = ■ = — - — { ) 

0 4H Q 4 3 x v V 0 ' 

where B Q is the radius of curvature of the streamlines 
and ~w 0 the z component of the velocity induced at the 

upper wing "by the lifting and trailing vortices of the low- 
er wing. The equation is based on the practically adequate 
assumption (see, for example, Fistolesi, reference 46) that 
a piano wing in a flow concave downward possesses the same 
lift as a wing convex downward possesses in an uncurved 
flow. Limiting ourselves to the linear members in tha dis- 
turbance velocities, we obtain 

t„ a w n 
0 4 ? a x 

It is of interest to compare the curvature effect obtained 

if in w n only the effect of the lifting vortex is taken 
u z 

into account and if also that of the trailing vortex is con- 
sidered. The result appears in all the computed examples 
which include the normal range of stagger-gap ratio a/h 
that the effect of the lifting vortex on the curvature of 
the streamlines is greater by far than that of the trailing 
vortex (fig.. 6). It can thus be understood why the previ- 
ously carried out computations of the curvature effect for 
two-dimensional flow where only the influence of the lift- 
ing vortex was taken into account gave relatively good re- 
sults. All the above considerations may correspondingly be 
applied to the lower wing (fig. 3). 

Before writing down the equations for computing the 
circulation, we shall introduce a few more brief notations: 
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for the upper wing, and 




for the lower wing. 

As may be easily seen f, , f„ , and f-, , f„ 

o "c u u . 

represent geometrical distances (see figs. 2 and 3), which 
may be read off from a drawing of the cellule and- thus the 

computation simplified. The terms f3 and f_ , respec- 
ts J u 

tively, were obtained by the combination of other expres- 
sions. 

Combining, we obtain from (1) and (2) after some trans- 
formations for the circulation distribution of the upper 
wing: 
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Vo 



0 



+ 1 



= a 0 - 



3 TT./ 

-1 



d 7 0 d TV 



d T) 0 > Tl 0 ' - \ 



(monoplane equation) 



u u 
+ 

8 TT L 



+1 I 

d r^ 1 

a 0 h ; 7 n 



<u 



— l 



(change in longitudinal 
velocity) 



+ i 



d Tl. 



+ a / 7,, ^— (downwash due to lifting vortex)- 



+ i 



- 1 



f 3 (a 2 + h 2 )- 2 f-i' 
3 o 



b.. D 0 \ 



4 > 



x o a o 



r<* ^u' 



(downwash due .to trailing vortex) 



3 h 

— - \ d t,.. 1 

3 f 5 J Tl 



(curvature effect) 



where 7 = T /bV . 

" • For the lower wing we obtain from equation (2): 
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7 Id 
'.xl ■ u, 



.1, 



C l u *T1 ^ 2 TT 



, 2 

b o r 



+i 

P .d .7 



u 



-l 



d V 'Ha 1 - % 



(monoplane equation) 



- a 



d r> 1 

a '0 



f x 3 

u 



(change in longitudinal 
velocity) 



0 f. 3 

■i ^u 



(downwash due to lifting vortex) 



/ 7e JT* + ' 

-i 2 u 



f 3u (a 2 + h 2 )- 2 (^ 0 .L^y 



f 3 f 4 

x u s u 

(downwash due to trailing vortex) 



i 



d Ti, 



+ "/ 7 °tf 

""- 1 



.? h | 



dT) 1 



3 5 ^> 'o 

x u x u 



(curvature effect) 



The sketches, (figs. 4 and 5) serve to clarify the 
signs of the angles of attack induced by the interference. 

It is particularly worthy of note in these equations 
that the parameters do not always enter linearly. The 

.term which characterizes the change in the longitudinal 
velocity is square with respect to the angle of attack (in 

: 7 Q there occurs a Q as a factor). The ' advantage s • of the 

linear theories, for example, the possibility of obtaining 
the total circulation through the superposition of the 
simple basic distributions are thus lost. 

. For our computation, we have only made, the assumptions 
which are made. for the individual wings of the Prandtl air- 
foil theory. For' the determination of the angles of attack 
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induced by interference, the cellule may be of arbitrary 
shape.. The two equations lead to the following conclu- 
sions as to the mutual interference of the wings of a bi- 
plane cellule: 

1. In general, the interference effect decreases 

with increasing stagger h; 

2. At the upper wing increases with increas- 

ing span of the lower wing; 

3. At the lower wing there similarly is a larger 

ai^ o corresponding to a larger span of the 

upper wing; 

4. If a > 0 (forward camber) almost all factors 

tend to a decrease in the effective angle of 
attack at the lower wing, whereas at the upper 
wing various effects may cancel each other. 
It is even possible with larger stagger to ob- 
tain an increase in the effective angle of at- 
tack for the upper wing. Forward stagger al- 
ways affects the lower wing more than the up- 
per. 



5. If a< 0 essentially the reverse of the above 
holds true. 



In order to obtain an idea as to the order of magni- 
tude of the various factors each is plotted for the upper 
wing of a recent biplane (example III) in figure 6, where 
there is also included a plot of the self-induced angle of 
attack . * 

The combined effect of all the factors results in 

downwash over the greatest portion of the span, so that 

CC; , according to equation (1), is mostly positive (in 
o u 

fig. 6 the po si t ive • direct i on is downward). The greatest 
effect is that shown by the tip vortices for which the 
part corresponding to the unstaggered cellule (term without 
a in (4)) is most prominent. The effect of the trailing 



*With regard to the peaks that occur in the cCj_ curve, 
it is to be observed tiiat these are due to a discontinuity 
in the twist of the upper surface. 
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vortices on the curvature of the streamlines is slight as 
was mentioned above, The- .curvature is mainly to be as- 
cribed to the effect of the lifting vortex. The lifting" 
vortex also exerts an effect through a change in the lon- 
gitudinal velocity and through up.wash although these efr 
focts are ' relatively small. From figure 6 the important 
fact is to be noted that, in a more accurate theory all 
terms in equations (4) and (5.) must be taken into account.* 

For determining the required 7 0 and y t we thus 

have the two integr odif f erential equations (4) and (5), 
whose particularly inconvenient property is that it is not 
possible to effect a separation into two equations for 7 n 
and 7 U . 

IV. METHOD OF SOLUTION 



A solution of these two. equations for 7 0 and 7 U 

in closed form is hardly worth attempting. *T.e shall thus 
try to determine the circulations from equations (4) and 
(5) by an approximation method where we may assume that 
the interference and self-induced angles of attack are 
small compared with the others. 

We first solve the equations for the monoplanes 
7 0 ^ 0 1 T +1 d y o 4 Tlo' 



J _o_ 

o.o 



= a. 



o 




'_ 1 d V V - 



and 



7 u *u = a _ J_ A ^u d V 

c i 11 t u U 2 tt ./ d V r u' " ^u 



-1 



♦That the sum of all the parts in figure 6 approximately 
agrees with the effect of the trailing vortices. is to be 
ascribed to the special shape of cellule III and cannot be 
assumed as true in general. 
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by the method of H. J/iulthopp (reference 47) and obtain 

from these the zero approximat i on . . 7 0 ^ and 7 U ^°\ 
"With these values, we next compute f rom e qua t i on s (4) and 

(5), a, (7 U (0) ) and a ± (7 0 (0) ) (the integrations are 
ou uo 

performed graphically) and substitute these values in the 
equations as a twist, that is, to be added to a Q and CC U , 

and from these new values and 7,, ^ ' are deter- 

mined, and thus the biplane interference effects accounted 
for to a first approximation. Similarly we obtain 7 n v2 ' 

and 7-u' 2 ', etc. In this manner a well converging- itera- 
tion process is made possible which gives us values for 7 Q 
and 7 U sufficiently accurate for practical purposes even 

after the first, but at most after the third approximation. 
Having determined with sufficient accuracy the values of 7, 
we may compute the induced angle of attack o.j_ , the total 
lift c*o ef f i ci en. t s of each wing' c a , .the induced drag coef- 
ficient c wi , the rolling moment coefficient c m q," and 
the coefficient for the induced pitching moment . c Iftgi for 

each wing and for the cellxile in the usual manner, as well 
as the remaining magnitudes depending on the lift distribu- 
tion'(see section II). 

As regards the practical carrying-out of the computa- 
tion the following may be said. The integrals of equations 
(4) and (5) are made up of a. few components that for the 
most part may be obtained graphically (see figs. 2 and 3). 
The determination of the approximations is particularly 
simple, since in the system of equations to be solved by the 
method of Muithopp, only the terms involving the twist por- 
tion vary. All the computations may be carried out by em- 
ploying technical aids. The time required is about a week 
for simple biplane arrangements up to three weeks for more 
complicated arrangements. Eiperienced computers can still 
further reduce the time considerably. 



IT.' GENERALIZATION' "OE' THE METHOD ' 

Up to the present we have assumed that each wing may 
be replaced by a lifting vortex" filament. It appears, how- 
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ever, that this assumption is not sufficient in some par- 
ticular cases. The first difficulty occurs in the deter- 
mination of 'the stagger a. We had assumed that the cen- 
ter of pressure at each section of each wing lies on 
straight lines and denotes the distances between these two 
lines "in the x direction by a. It is evident that in 
the neighborhood of the zero lift of a wing this method 
leads to contradictions since the center of pressure then 
travels off to infinity. Furthermore', very close to a 
wing the velocity field due to the substitution vortex 
does not agree sufficiently with that of a wing of finite 
chord, thickness, and camber. We are thus constrained to 
improve the theory given above for the range near the zero 
lift and the case of wings of small gap. A comparison of 
the computations and measurements given below will show, 
however, that the simple theory with a single vortex ele- 
ment is sufficiently accurate in most cases if the gap- 
chord ratio is equal to or greater than 1. To present the 
relations correctly it would be necessary to replace each 
wing by a lifting surface. A step in this direction is 
the following method, which already gives sufficiently good 
results. ' * 

We replace each wing by two straight-line vortices at 
right -angles to . the flow and having strengths T 1 and T 2 , 
the position of the vortices being assumed at-l/4 and l/z 
chord, respectively, from the leading edge of the wings.* 
This method has the following advantages. In the first 
place, by using two vortices the velocity field is more 
accurately represented. In the next place it is now pos- 
sible at each position along the span to have the center 
of pressure of each of the substitution vortices agree 
with the actual center of pressure at each wing section by 
suitable distribution of the circulation of the two vortex 
elements. If only a single lifting vortex element is as- 
sumed, it is impossible to take into account deviations of 
the centers of pressure from the straight lines. Finally, 
on the assumption of two vortex elements, it is possible, 
in addition to effects on the angles of- attack, also to 
consider more accurately the curvature of the streamlines 
due to the interference. Since the first vortex element 
is assumed at l/4 chord and the second at l/2 chord from 
the leading edge, the forward element, according to the 



*Iiore general singularity substitutions for the two-dimen- 
sional biplane problem will be found, for example, in B. 
Pistolesi (reference 30). 



IB 
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-theory, will take up the effects of the angle; of attack and 
the rear element, those of the streamline curvature. Prac- 
.tically, however, these effects are not so easily separated. 

We must first set up" a theory of the monoplane wing 
which is replaced by two lifting vortices-. It is to be 
noted, however , that this theory is of use only in regard to 
our final object since in the case of the monoplane n.lone 
the self-induced curvature of the streamlines is so small 
that it may be neglected. In addition to the' -angle of at- 
tack equation 



ct. 



ct„ - 
g 



a. 



(6) 



we have the relation 



= f - f 



(?) 



for the camber f denoting in general the rise in camber 
of the section. In these equations we must express Ct e , 
f e , and ctj , fj through the circulations of the two 
lifting vortices. To determine ct e and f - e , we consider 
the relations at a wing cross section. For the lift and 
moment coefficient the following^ relations hold with good 
approximation 



with 



- c 1 

- c a 



c ' ' 
a 



Ct e + 



°m = °ra l (cte + 



/ d , c a \ 
~ V d a J 



* f e 



A 



) 



> (3) 



= ( 



d c r 



v d a 



A — > 



Me must further express c a and Cjjj in terms .of the cir- 
culations r x and T a of the lifting vortex elements in 
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order that we may have the relation between cc e » fe, and 
1^ , P 2 . A relation between these magnitudes is given by 
the requirement that the strength of' the vortices Y x ' and 
T 2 be such that the total lift and total moment at each 
wing' section is not changed by the double vortex substitu- 
tion.* 

r x + r a = X V t (9) 
4 ^ + I P 2 = T 7 t (10) 



With the aid of these relations and equations (8), there 
is obtained 

a e = 2 -±- - - Z— (11) 

V 2 V 
V t V t 

* 2 

2 c a ' 4 c > 

v t -2- v t 



In equations (5) and (?) there are still to be determined 
and i^. Tor this we required the sum of all induced 

downwash velocities w at the midpoint (x = - t/4) and 
we obtain 

, •■ f i t 1 3(w/V) 
. a± = w/7 and 



t 8 R .83 (x/t) 



*The horizontal components of the disturbance velocity, 
which produces F 1 at the position of F 2 and F a at 

the position of F x is assumed small in comparison with 
V and its effect is therefore neglected. 
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where E is the radius of curvature. at the center of the 
wing section through self-induction of the curved stream- 
line. The value of w is determined by & method of sec- 
tion III, which need, hot "be gone into here,. We obtain 
finally from equations (6) and. (7) with (11) and (12) , set- 
ting 7 = r/b 7: 



2 b 



i y x + 2 7 a 

4 



G_ ' t 



2 b 



2 TT 



+ 1 



d -n' 



i + 



2 b 




t y 



— J + On - -n») L_ 

2 b/ _ rjt 



.+ 1 



d 7 a an» 



d iy ri - iy 



(13) 



7j + 7 2 1 y + 2 7 £ 



+ - 



c a i t 4 c m « t 



2 b 



2 b 



2f g b 



+i 



Hi 

d Tl 



+ 



4 TT b 



+1 



d T.« 



-1 



fc Y + («n - TV 




a ti» (14) 



2 b' 



There are thus set up the monoplane equations based upon 
the assumption of two vortex elements for each wing. 
These equations may be solved by iteration by determining, 
for example, as a rough approximation the sum y + 7 S 
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according to the usual airfoil theory, separating into 
Tfj -and 7 2 according to the two-dimensional theory (equa- 
tions (9) and' (10)) and then using these as initial values 
for the iteration process, substituting in the right sides 
of equatidns (13) and (14). 

Through the biplane interference effects, the longi- 
tudinal velocity, the angle of attack, and the curvature 
of the flow are changed by amounts AV 2 . A oc 2 , Af 2 , and 
from these there is obtained a change in and "r g by 

amounts Ar, and AP- given by* 



A V a .3 r 

Ar = r l — 3 - + — 

1 1 V 9a 



i r 9 r i 

i- A a 2 +1 - A f 

L 3 f 



Ar ; 



A V, 



+ 



s r 3 

L 3 a 



Aa 



2 ! + 

J 



3 f 



A f. 



These equations may also be put in the form** 



A Vg „ 

Mi = Vx ~ + c ai - -1- A a 2 + I 



A7 a = 7 P 



x 2 o 



a v 3 r t 



' 2 2 V ' L ~ 2 2 b 



j ^' A f 3 
L b 



Act J + 



cap. ' A f £ 



> (15) 



In these equations we determine c ai ' . and c a. Q * *>Y plot- 
ting first c a (a =o ) for the entire wing and separating in- 



*Since according .to the ■ theory c m ' = c a l /4 
the square brackets drop out because 



3 r 



the terms in 

3 r. 



1=0 and 



3 f 3 cc 

= 0.. . Actually this assumption holds generally only within 
a restricted range of angle of attack. It is necessary to 
carry these' terms along in some cases, therefore. 

**It is here assumed as above' that the familiar relation 
c a = c a ' (cc e + 2 f e /t) holds for each wing, see formula (g) 
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to c a and c a according to equations (9) and (10) ; • 

The slopes of. these curves give us the required values. 
It will then appear whether or' not c a = 0, i.e., wheth- 
er or not the separation was necessary. If we not substi- 
tute in these equations for the disturbance values, the 
expressions obtained in section III and also in equations 
(13) and (14) 7 \ + A 7 x and 7 2 + A 7 2 in place of Y x 
and 7 2 , respectively, we have the new integrodif f erential 
equations, which must be satisfied by 7 X and 7 2 . There 
are obtained altogether four equations for the four un- 
knowns 7 Q , 7 D , 7 U , and 7 U . 

° i °a i s 

VI. RESULTS OF THE COMPUTATION 



We shall now make use of the above method to compute 
the lift distributions of several types of airplanes. Of 
the five models computed (fig. 7), measurements on only 
the first three were available for comparison. Models IV 
and V were freely chosen, the former representing the con- 
ventional biplane cellule (two equal rectangular wings 
with rounded-off corners of aspect ratio 6 and gap-chord 
ratio 1). The computation was performed for zero and mod- 
erate forward stagger of the upper wing. Model V repre- 
sented the limiting case of the sesquiplane with ellipti- 
cal wings. In presenting the results, there was first 
computed the total lift coefficient c a from the formula 



c a = A / 7 d r, ■ 

•/., 

for each wing and for the biplane arrangement and the val- 
ues plotted against a. The angles of attack induced by 
the mutual interference and the spanwise lift distributions 
are given in section VII below where a comparison is made 
of the various models. Unfortunately, only in the case of 
model I was it possible to compare the lift distributions 
with the results of measurements. 

Example ;1; The computation was first carried out for 
the cellule of a rather old type airplane* (model I), the 



♦Shape of cellule and other data as in the case of the re- 
maining examples are given in figure 7. 
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American single-seat' pursuit airplane' PW-9 for which' de- 
tailed measurement s were "available (references 38 and 39). 
Thi s.. example was chosen "because pressure distribution; 
measurements on each wi.ng and on the biplane arrangement 
were . avai lable so that the computation could be compared 
with the results of tests. In figures 8-and 9 the comput- 
ed and measured distribution of 7 = c a t/2 b was plotted 
against the span for each wing for various angles of at- 
tack. It is seen that there is good agreement between the 
two. On the upper wing the first approximation is already 
sufficiently accurate while for the lower wing, for which 
the deviation from the undisturbed wing is considerably 
greater, a second approximation was necessary. A third . 
approximation would only slightly affect the results. The 
total lift is plotted for each wing as a function of 

the angle of attack in figures 10 and 11. The indicated 
test results of Fairbanks and Loeser (differing somewhat 
from each other) show that the computation lies entirely 
within the range of measuring accuracy. 

Example 2: iiodel II (fig. 7) is a biplane cellule 
dating from the G-ottingen measurements at the time of the 
World War (reference 28). We were led to measurements on 
this model as a result of a still unpublished work by E.. 
Fuchs, where for this particular example no satisfactory 
agreement in the valxxes between the computation and 

measurements could be attained. This model has the pecul- 
iar property of having the lower wing set at 4° less than 
the upper wing. Moreover the gap-chord ratio is quite 
small (h/t = 2/3). With regard to the measured value's 
there is to be observed the following. In the case of the 
old Gottingen measurements, the stagger could not accurate* 
ly be held at the constant value 0.5 t on which the compu- 
tation of Fachs and also our own computation was based. 
Thus the measured values for the lower wing are derived' 
from the value a = 0.5 t and those of the upper wing from 
the value a= 0.547 t. These measurements were then re- 
computed in terms of a =0.5 t (in fig. 12 both values are 
plotted) so that some uncertainty is attached to the points 
indicat ed . 

In carrying out the computation, it appeared that' the 
method involving a single lifting vortex filament for each 
wing was not -adequate, the main reason for which was to be 
found in -the small gap of the biplane. In addition, a re- 
finement in the computation was necessary for the cellule 
a.ngle-of-attack range in the neighborhood of zero degrees. 



24 



N.A.C.A. Technical . Memorandum So.. 889 



In this range the lower wing has practically.no lift While 
the measurements on the upper wing, which. was still 4 9 re- 
moved from the zero lift'angle showed an effect "by the 
presence of the lower wing (figs. 12 and 13). We therefore 
made use of the method'given in section V above. We first 
replaced the lower wing in the angle-of-at tack range ■ in ' 
question by two lifting vortex elements while' still using 
a single vortex element for the upper wing.' The result •■. 
was that also at the upper wing at a '= 0° there was. an 
increase in lift and at the lower wing a decrease in lift 
as compared with the values obtained on the assumption of 
a single vortex for each wing. The agreement with the 
measured values was still not quite sufficient. We there- 
fore attempted to improve the computation by the assumption 
of two lifting vortex filaments for each wing. This gave 
an increase in lift for the upper wing over the entire 
angle-of -attack range, so that good agreement was obtained 
with the measurements. Sufficiently good agreement was in 
this manner obtained also for the lower wing. 

Example 3: The c ellule of a . modern type airplane* 
(fig. 7, model III) served as the next example. Figures 
14 and 15 again show the measured. and computed decrease in 
the total lift coefficient for the upper and lower wing 
and it is seen that the computation sufficiently accurate- 
ly represents the conditions. The deviation from the unr-.. 
disturbed wing is in this case not as large as in the pre- 
vious examples. A more exact comparison of. the three, 
cases will'be given below (section VII). In the neighbor- 
hood of zero angle of at tack . of .. the wing, the method given 
in section V was used, each wing being, replaced by two 
lifting lines. For this- example, we also computed- the lift 
distribution for each wing for several deflections of the 
ailerons, which were located on the upper wing. Figure 16 
shows the externally induced angle of attack . a i uo at the 

lower wing' for an aileron deflection j3q = ±1 (pq = +1 

on the' right wing and -1 .on the left). The lift distri- 
bution of the upper wing is shown in figure 17 and that of 
the lower wing 'in figure 18 for various aileron deflections 
ang angles of attack. The effective rolling-moment coef- 
ficient c • for the entire cellule would be, without in- 
m q • 

terference, +0.582 for |3 ^ = ±10° and + 1.166 for. = 

±20°.. On account of the modified lift distribution of the 
lower wing, an effective rolling moment in the opposite di- 
rection is produced so that the total rolling-moment coeffi- 

*The test results were kindly made available to us by the 
manufacturer s . 
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cient of the "biplane cellule' is only +0.564 for |3q Q = 

±10° and +1.13 for fiq^ = ±20°. (The angle of attack 

enters practically linearly and therefore has only a 
slight effect on c^, ) Measurements of this kind that 
could be used for comparison were not available* The re- 
lations' given are naturall-y valid for conditions "below the 
flow separation point.'- The effect of the change in the 
lift' di stribution- as compared with that of the monoplane 
on the separation will be discussed later (section VII). 

Example 4: The computation for- model 17 (fig. 7) was 
carried out once for the unstaggered biplane (a = 0; IVa) 
and again for the staggered biplane (a/h = +0.45; IVb) by 
the method of section IV. Since the upper and lower wings 
are equal, the computation, was very simple as many expres- 
sions could be determined with the two wings. in common. 
Figures 19 and 20 show the total lift coefficients as a 
function of the angle of attack for each wing and bring out 
the effect of stagger. The result -showed that the upper- ' 
wing of the staggered cellule experiences hardly any de- 
crease in lift, as compared with the unaffected wing, -rhilft 
the decrease is stronger for the lower wing, so that as a 
whole the cellule experiences a decrease in lift. Figures 
19 and 20 further show that the third approximation is al- 
so sufficient for this case where the mutual interference- 
is relatively large. In figure 21 are shown the polars* 
for the upper wing, lower wing, and cellule, which, in case 
of equal wi ngs must approximately agree. This is at the 
same time a check on the accuracy of the measurements and' 
shows that in spite of the many graphical integrations the 
computation was carried out with sufficient accuracy. 

Example. 5: To complete the series of examples, there 
was computed as a final example a cellule for which the span 
of the lower wing was half as large as that of the : upper 
wing, i.e., a sesquiplane (model V, fig. 7). The wings of 
the freely chosen cellule are elliptical. Further, the gap 
between the wings was chosen relatively small (h/t 0 = 0.75). 
In the rather large resulting interference effect, which 
changed the elliptic lift distribution considerably, it was 
necessary to -carry out the computation for each wing as far 
as the third approximation. Figures 22 and 23 show, however, 
that a further refinement is not necessary. , 

*c^T was determined from the li f t di str ibuti ons from ■ c^ 

= A / 7) d r\. 

-i 
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VII. CQMPARI SON . OF THE .VARIOUS EXAMPLES 

We shall now compare the results for the various mod- 
els with regard to several aerodynamical characteristics. 
For this purpose there -was plott.ed on figure 24 the inter- 
ference angle of attack at the upper and lower wing of each 
'model for the angle; of attack 5° (above zero angle of at- 
tack) of the cellule so that the mutual interference ef- 
fects could be compared in. each case. Further, in figures 

.25 to 29 are plotted the ratio — , the chord t, and 

the ratio-' of the local c a values to the total c a value 

of the wing under consideration. Figures 25 to 29 there-, 
fore give us information as to the position along the span 
at which there is overload- and underload, respectively, 
and thus show where the flow will first separate in the 
neighborhood of c a . From the standpoint of rolling 

stability, it is de sirable" that, the middle portions of the 
spans have a larger iocal c a than the outside portions, 
since in that Case the separation will start in the neigh- 
borhood of the fuselage. and spread out toward the wing tips. 
Otherwise, separation starts at the wing tips, a result 
which may easily lead to loss in lateral stability. It is 
readily seen that an upwash induced at the wing tips due to 
interference opposes a favorable lift distribution without 
interference effect. The effect will be unfavorable in the 
neighborhood of ^a max part icularly- with regard to the ef- 
fectiveness of the ailerons which in general operate in 
this range. Such an upwash region is found in the case of 
the -upper wing of model I, figure 24. The corresponding 
range of overload in figure 25 is shifted by the presence 
of the .-lower wing farther • toward the wing tip. Since also 
the lower. wing of this model ' show s the . same effect -to -a 
small extent, it is to be supposed that the entire cellule 
will have at large angles of attack unfavorable character- 
istics -with regard to the lateral stability. 

Since the upwash region discussed above is to be as- 
cribed to the effect of the trailing vortices, model II 
(fig. 26) and model. I?, on account of b Q = "b , cannot 

show 'this effect ,(fig. 24). The upwash at the upper wing 
of model II is due to the lifting vortex. The lower wing, 
however, receives a downwash of such magnitude that in 
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order to attain again the lift that it would have in the ,. 
undisturbed flow, it would have. ,to,be set over a large . 
part of the span, at an angle 5° to 6° higher. Model II, 
however, does not show the unfavorable lateral stability 
characteristics of model I, a fact that is to be ascribed 
to the rectangular shapes and the equal spans of the two 
wings. 

The most favorable characteristics in every respect 
are shown by model III. The mutual interference as com- 
pared with that of the other models is smallest (fig. 24). 
In spite of the fact that the lower wing of model III as 
well as that of model I has a smaller span than the upper 
wing (bu/bo = 0.763. for model III and 0.708 for model I), 
the upwash at the tips of the upper wing of model III has 
to some extent been avoided by suitably shaping the wings 
so that the good load relations (fig. 27) are not unfavor- 
ably changed to any large extent and it may be stated that 
this model shows good characteristics with regard to lat- 
eral stability. 

With regard to model IV (fig. 28), it is worth noting 
that the overloading due to the interference effect is 
somewhat equalized as appears particularly in the lower 
wing of the staggered cellule IVb where c & t /c a and t 
almost agree. The upper wing of model IVb shows almost 
no effect (fig. 24). Otherwise, model IV shows character- 
istics similar to model II. The downwash at the upper wing 
of model IV, however, due to the larger gap h of the 
wings (fig. 7) is not so large (3° to 4° over a large part 
of the span) . 

The particular reason for carrying out the computation 
on model V was that a larger upwash region was expected at 
the upper wing. This expectation was confirmed by the co-m- 
putation (fig. 24). The effect, however, is not so large 
as not to be avoidable by a correspondingly moderately 
large twist. There is to be sure such a large change- in 
the elliptic lift distribution (fig. 29) that the separa- 
tion begins at t.he wing tips so that this model does .not 
appear to have safe lateral stability characteristics. For 
biplane arrangements such as those of these examples, it 
may be said that at least for the upper wing the rectangu- 
lar plan form has advantages aa compared with others.. 

It may be pointed out in concluding that in addition 
to the lift distribution, it is also possible to compute 
the induced drag as has been done for one case (fig. 21). 
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A comparison with measurements gives. the noteworthy result 
tha t (for example, for model -III) -the profile drag c w - 
c w ^ (where .TZfi is .taken froia measurement a), was quite con- 
siderably increased at the' lower, wing through the effect 
of the upper wing and shoWed a dependence on the angle, of 
attack (fig. 30). To explain this it' would be necessary 
to investigate the pressure and velocity relations in the 
boundary layer of the lower wing, which, however, would lead 
us beyond the confines of this report. 

VIII. SUMMAEY 



After a brief review of the work that has so far been 
done on the biplane problem, there are set up the general 
equations for the circulation distribution along the span 
of the wings of a biplane cellule and a' method of solution 
indicated. The assumption is first made that each wing may 
be replaced by a vortex element in the senso of the Prandtl 
airfoil theory. In the examples computed it appeared that 
this assumption was' sufficient for the usual biplane ar- 
rangements and the results of the computation show good 
agreement with the measurements. A refinement of the the- 
ory is necessary, however, in the neighborhood of zero an- 
gle of attack and for the case of very small gap-chord ra- 
tio. By the assumption of a double vortex system at each 
win;;, better agreement is obtained and the curvature of 
the streamlines is likewise better taken into account. 
With this' assumption, there is first set up a theory of the 
monoplane and this theory is then applied to the' biplane, 
fith the aid of five examples, which cover the range of us- 
ual biplane arrangement s, '" i t is shown how -large is the ad- 
ditional angle' of attack induced at each wing by the pres- 
ence of the other. • The examples are also compared with 
regard to stability about the longitudinal axis in the 
neighborhood of' ^a" max « "' " ' 

The author i s very grateful for the many valuable sug- 
gestions of Miss I : . Lotz- during the. conduct of this inves- 
' ti gat ion .' • • ' •• 
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Figi. 1.2,3,4.5,6 




Figure 1.- The vortex system 

used as a basis 
in the investigation. 
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Figure 2.- Sketch for obtaining the 

relations at the upper wing. 




Figure 3.- Sketch for 
obtaining 
the relations at the 
lower wing. 
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Effect of change in longitudinal 
velocity by bound vortex. 
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Figure 4.- Sketch jig^e 5.- Sketch 

showing showing 
various conditions various conditions Figure 6„- Distribution of oV 
at upper wing. at lower wing. for model III. 



a~ Total angle of attack induced 
by the mutual interference 
effect . 
b- Part due to tip vortices 

without stagger, 
c- Part due to the stagger of the 

tip vortices. 
d» Part due to the curvature pro- 
duced by the lifting vortices, 
e- Part due to the curvature by 
^ the trailing vortices. 
*V f— Part due to the lifting vortices 
(vertical component), 
g- Part corresponding to the 
change in the longitudinal 
velocity. 
Below; self induced angle of 
attack for comparison. 
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ZigB. 7,8,9 
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Figure 7.- Computed examples. Shapes 

of cellules, profile 
sections, and wing setting referred 
to chord of middle section in upper wing. 




a- Wing without interference, b- Wing with mutual interference 
effect, first approximation. c-Wing with interference, second 
approximation.»-Wing without interference measurements results, 
O - Wing with interference effect, measurement results. 



Figure 8.- Distribution of 

circulation over 
the half span* 



Figure 9.- Circulation distribution 
oyer the half-span. 
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Figs. 10,11,12,13,14 
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a- Wing without mutual 
Interference, b- Wing 
with interference, first 
approximation, c- Wing 
with interference, second 
approximation. • Wing 
without interference, 
measurements of Fairbanks. 

X Wing without interfer- 
ence, measurements of 
Looser. o Wing with inter- 
ference, measurements of 
Fairbanks. » Wing with 
interference, measurements 
of Loeser. 



e . .t « w 

Figure 10. -Lift curves. Figure ll.-Lift curves. 
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a- Wing without interference, computed, b- Wing 
with interference, computation with one vortex 
at upper wing, two in lower, c- Wing with inter- 
ference, computation with two vortices in each 
wing. • Wing without interference, measurement. 

9 Wing with interference , uncorrected measure- 
ment ( a =0,54? t ). o Wing with interference, 
corrected measurement ( a=0.5 t ). 

Figure 12.- Figure 13.- Lift curves. 

Lift curves. ( The break in 

curve a. is only slightly 
noticeable in curves b 
and c) . 
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a- Wing without inter- 
ference, b- Wing with 
interference. • Wing 
without interference, 
measurement, o Wing 
with interference, 
measurement. 

Figure 14.- Lift curves. 
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Pigs. 15,16,17,18,19 
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a- Wing without interference, 
b- Wing with interference, c- 
Wing with interference, second 
approximation. • Wing without 
interference, measurement. O 
Wing with interference, 
measurement. 

Figure 15.- Lift curves. 
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Figure 16.- Magnitude of interference 

angle of attack at lower 
wing for deflection of aileron pn x 1, 
for various angles of attack of cellule. 




a- Without interference by the lower 
wing for Pq 0 = 10° and 200. b- with 
interference by lower wing, /3q^= 10o. 
c- with interference by lower wing. 
13 Qo- 20°. 

Figure 17.- Circulation distribution 

for aileron deflections 
at upper wing for various angles of 
attack. 



a- Wing without interference; 
b- wing with interference, un- 
staggered, first approximation; 
c- wing with interference, un- 
8taggered, second approximation; 
d- wing with interference, un- 
staggered, third approximation; 
e- wing with interference, 
staggered. 

Figure 19.- Lift curves. 




a- Without interference; b~ with inter- 
ference, |3q = 0°; c- with interference, 
PQp'lO 0 ; cl- with interference tfcfo- 20°. 

Figure 18.- Circulation distribution 

for the lower wing for de- 
flection of ailerons at upper wing. 
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ft- Wing without interference; 
b- wing with interference, un- 
staggered, first approximation; 
c- wing with interference, un- 
staggered, second approximation; 
d- wing with interference, un- 
staggered, third approximation; 
e- wing with interference, 
staggered. 

Figure 20.- Lift curves. 
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Tigs. 20,21,22,23 

J$fSL 




OK Opt 

a- Upper and lower wing without 
interference; b~ upper wing with 
interference, unstaggered. c- 
upper wing with interference, 
staggered; d- lower wing with 
interference, unstaggered; e- 
lower wing with interference, 
staggered; f- cellule wings 
alone; g- cellule unstaggered; 
h- cellule staggered. 

Figure 21.- Polar diagram. 




a- wings without interference; b- wings with interference, first 
approximation; c- wings with interference, second approximation; 
d- wings with interference, third approximation. 



Pigure 22.- Lift curves. 



Pigure 23.- Lift curves. 
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Tigs. 24,25,26 
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Figure 24.- Comparison of angle 

of attack induced 
Dy biplane interference at both 
wings ( a. = 5° from sero angle 
of attack). 



a- wings without interference; 
b- wings with interference. 



Figure 25.- Lift distributions, j 
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a- wings without interference; 
b- wings with interference. 

Figure 26.- Lift distributions. 

The single curve for 
the upper wing holds for the wing 
with, as well as without inter- 
ference. 
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figs. 27,28,29,30 




a- wings without interference; 
b- wings with interference. 

Figure 27.- Lift distributions. 
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a- wings without interference; 
b— wings with interference, un- 
staggered; c- wings with inter- 
ference, staggered. 

Figure 28.- Lift distributions. 
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a- wings without interference; 
b- wings with interference. 

Figure 29.. Lift distributions. 
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a- wings without interference; 
b~ wings with interference. 



Figure 30.- Profile drag of wings. 



